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Abstract: If the cold dark matter consists of weakly interacting massive particles 



(WIMPs), anticipated measurements of the WIMP properties at the Large Hadron 
Collider (LHC) and the International Linear Collider (ILC) will provide an unprece- 
dented experimental probe of cosmology at temperatures of order 1 GeV. It is worth 
emphasizing that the expected outcome of these tests may or may not be consis- 
tent with the picture of standard cosmology. For example, in kination-dominated 
quintessence models of dark energy, the dark matter relic abundance can be signifi- 
cantly enhanced compared to that obtained from freeze out in a radiation-dominated 
universe. Collider measurements then will simultaneously probe both dark matter 
and dark energy. In this article, we investigate the precision to which the LHC and 
ILC can determine the dark matter and dark energy parameters under those circum- 
stances. We use an illustrative set of four benchmark points in minimal supergravity 
in analogy with the four LCC benchmark points. The precision achievable together 
at the LHC and ILC is sufficient to discover kination-dominated quintessence, under 
the assumption that the WIMPs are the only dark matter component. The LHC 
and ILC can thus play important roles as alternative probes of both dark matter and 
dark energy. 
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1. Introduction 

Current cosmological data [1-4] indicate that the energy density of the universe 
today is dominated by degrees of freedom beyond those of the Standard Model (SM) 
of particle physics, with approximately 73% given by dark energy, a form of energy 
characterized by negative pressure, and approximately 23% given by nonbaryonic 
dark matter. Given that future colliders such as the Large Hadron Collider (LHC) 
and International Linear Collider (ILC) are designed to probe physics beyond the 
SM, it is important to investigate whether and to what extent these experiments can 
probe these elusive forms of matter and energy. 

For dark matter, it is well known that both direct and indirect collider physics 
connections can be established, depending on the properties of the dark matter can- 
didate. One particularly well-motivated class of models is that in which the dark 
matter is a neutral weakly interacting massive particle (WIMP). Such particles are 
in thermal equilibrium in the early universe, and the standard freeze out calculation 
predicts a value for their relic density which is in the right ballpark as required by 
cosmology. WIMPs are naturally present in models which attempt to solve the gauge 
hierarchy problem, since such models generically introduce new weakly interacting 
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particles. Prototype scenarios include models with softly broken N — 1 supersymme- 
try (SUSY), in which there is a conserved discrete symmetry, R-parity. In such mod- 
els the lightest superpartner (LSP), typically a neutralino, is a viable cold dark mat- 
ter candidate. The connections between astroparticle and collider physics have been 
extensively explored in the literature on supersymmetric dark matter [5-22]. More re- 
cently, other WIMP candidates have emerged in new physics scenarios which include 
a conserved discrete symmetry, such as models with flat [23,24] or warped [25,26] 
extra dimensions, and Little Higgs models [27-31]. 

It is much more difficult to establish a direct link to collider physics for the case 
of dynamical dark energy, since it is typically far more weakly coupled to the SM 
(for an example of an intriguing exception, see [32]). Although this energy density 
can be the cosmo logical constant, in which case cosmology has given us an invalu- 
able clue to the cosmological constant problem, it is plausible that the dark energy 
can be associated with an effective scalar field degree of freedom, commonly called a 
quintessence field [33-37]. Collider experiments are unlikely to probe the unknown 
quintessence field directly, given its typical range of masses and couplings. Connec- 
tions between collider physics and dark energy in the form of quintessence are thus 
necessarily indirect, and only a few are known. One connection is in the context 
of models with low energy supersymmetry. Since four-dimensional de Sitter space 
does not admit unbroken N — 1 supersymmetry [38], the discovery of supersymmetry 
through measurements of spins and couplings, and hints for its spontaneous breaking 
mechanism through measurements of mass patterns, will give us clues to the cosmo- 
logical constant problem and possibly its associated quintessence dynamics. Other 
indirect connections with collider physics result because quintessence dynamics can 
leave an imprint on the cosmological history, which in turn determines an observable 
quantity that depends partly on parameters deduced from collider measurements. 

One intriguing collider physics connection can arise through the gravitational 
interactions of the dark matter with quintessence. For thermal relics such as LSP 
dark matter, freeze out occurs when the expansion rate becomes larger than the 
interaction rate for reactions which change the LSP number density. In the standard 
cosmological scenario, freeze out represents the gravitational interactions between the 
cold dark matter and the relativistic gas of particles dominating the energy density 
of the universe at the time of freeze out. However, if the coherent field energy density 
of the quintessence field dominates during freeze out, freeze out now represents the 
gravitational interactions between the dark matter and the dark energy field degrees 
of freedom. In such cases, the relic density of the WIMP can be very strongly affected. 
More explicitly, the WIMP energy density today can be written as 




(1.1) 



2 



where Hp is the Hubble expansion rate at freeze out, m x is the mass of the WIMP, 
(<Jav) is the thermally-averaged annihilation cross-section, xp = T/m x ~ 1/20 for 
electroweak scale cross-sections (with only logarithmic dependences on Hp, (uav), 
and m x ), and T today is the temperature of the cosmic microwave background (CMB) 
today Eq. ( Jl.lp demonstrates that if Hp is increased due to the quintessence energy 
density, the WIMP relic density can be enhanced for the same microphysics con- 
tained in (o~av). Note that Eq. (|1 . 1|) also makes manifest the well-known property in 
standard cosmology that Q x h 2 becomes approximately independent of the mass m x 
for a fixed (cav), since in this case Hp ~ x F m x /M p i. 

A challenge for such alternate scenarios is that standard big bang nucleosynthesis 
(BBN) constrains any extra contribution to the relativistic energy density. The 
quintessence energy density therefore has to be large during the time of WIMP 
freeze out (T ~ 0(1) GeV) and then dilute more quickly than radiation to become 
effectively irrelevant by the time that BBN commences (T ~ 10~ 3 GeV). As pointed 
out by Salati [39], this can occur for scenarios in which the universe is driven by the 
quintessence kinetic energy (the kination-dominated period) before BBN. (Related 
scenarios were also suggested before by [40,41].) Long after BBN, the universe can 
enter a quintessence potential energy dominated regime. Such kination dominated 
freeze out scenarios are then consistent with standard cosmology and predict that 
the standard relic abundance computed from the parameters extracted from the 
next generation of colliders will be mismatched from the relic abundance deduced 
from observational cosmology. Indeed, the kination hypothesis is only one of many 
possibilities we must turn to if future experiments find a mismatch between the 
standard computations of the relic abundance and cosmological observations. 

The implications of kination domination for LSP dark matter have been explored 
previously [42-48]. We extend this work by analyzing the precision to which the LHC 
and ILC can simultaneously probe dark matter and dark energy. While colliders 
have traditionally been viewed as tools for discriminating among different particle 
physics models, we shall demonstrate that they are also capable of achieving sufficient 
precision on the effective field theory parameters to discern the difference between 
different cosmological models, e.g. the standard scenario and the kination domination 
scenario. The anticipated collider data will allow us to probe in the laboratory a 
new era in cosmological history, which may lead to significant deviations from the 
standard cosmological picture [49-52]. 

The order of presentation will be as follows. We begin by providing an overview 
of the issue of experimental tests of the WIMP hypothesis in Section ||. In Section [3], 
we review the effect of kination-dominated quintessence on the calculation of the 
WIMP thermal relic abundance, and provide a simple map between the kination 
scenario and the standard cosmological scenario. In Section |4] we then investigate 
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the capabilities of the LHC and ILC in pinning down both the dark matter and 
quintessence parameters. Following recent ILC studies, we use study points defined 
within the "mSUGRA" or "cMSSM" supersymmetric framework, where the super- 
partner masses are assumed to unify at the grand unification scale. Our study points 
are analogous to the set of four LCC study points (see e.g. [18]) which were chosen to 
represent the four "good" dark matter regions of the mSUGRA model. We reserve 
Section || for our discussion and conclusions. 

2. Colliders as dark matter and dark energy probes 

We begin by discussing in general terms how high energy colliders can test the 
WIMP hypothesis. The physics process which determines the present dark matter 
relic anbundance, is the dark matter annihilation into all possible final states with SM 
particles. Once we know the total annihilation cross-section a a of any given WIMP 
candidate, we can straightforwardly predict its relic abundance. Therefore, we can 
test the WIMP hypothesis by studying the properties of the WIMP candidates, and 
consequently, trying to constrain their annihilation cross-section a a- 

This is clearly a challenging exercise since colliders are not recreating the process 
of dark matter annihilation per se. Instead, given sufficient energy, colliders would 
produce WIMP dark matter particles, either directly or indirectly (in the decays of 
other particles). If those signals can be identified over the SM backgrounds, they 
can be studied in order to determine the properties of the dark matter particle, most 
notably, its mass and interaction strength to the different SM particles. It is in this 
sense that colliders are helpful in probing the dark matter. 

However, it is clear that the extrapolation from collider data to a prediction 
about the dark matter annihilation cross-section a a is rather challenging. The main 
problem is that the relic abundance is determined by a multitude of final states, 
while at colliders we typically observe the dark matter particles only in a limited 
number of exclusive final states. It is unrealistic to expect that colliders will be able 
to measure every single dark matter coupling, and the best one could hope for is that 
the colliders will be able to measure the largest (and therefore most relevant) cou- 
plings with some precision, while placing limits on the remaining (hopefully smaller 
and less relevant) couplings. Obviously, the way this is done is by either observing 
a dark matter signal in a specific channel, which would provide a measurement of 
the corresponding coupling, or failing to observe a signal in a specific channel which 
would provide an upper limit on the corresponding coupling. Then, by adding the 
results for all possible annihilation final states, one would obtain both a lower limit 
a™ and an upper limit a^ ax on the total annihilation cross-section a a- Assuming 
standard cosmology, these limits would correspondingly translate into an upper and 
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lower limit on the WIMP relic density Q x . From the preceding discussion it is clear 
that <7™ n (and the corresponding upper limit on Q x ) is derived from collider informa- 
tion from channels which have been observed while a™ ax (and the lower limit on Q x ) is 
derived from collider information from channels which have not been observed. It is 
therefore also clear that the collider information about <r™ n is more robust than the 
information regarding a™ ax . First, the presence of additional annihilation channels 
into non-SM final states would invalidate a bound on a™ ax which had been derived 
strictly within the SM. Second, arriving at a meaningful bound on a™ ax is experimen- 
tally quite challenging: the unobserved channels typically outnumber the observed 
ones, and also, the unobserved channels tend to be plagued with larger backgrounds 
and, correspondingly, larger uncertainties - after all, the large backgrounds are often 
the reason those channels have not been observed in the first place. 

With all this in mind, let us now turn our attention to some specific scenarios. 
For any given WIMP dark matter candidate cosmology data provides a preferred 
value a e ^ p for its total annihilation cross- sect ion. This value depends mostly on 
whether the dark matter particle is an s-wave or p-wave annihilator, and is almost 
insensitive to the dark matter particle spin [53] . The test of dark matter at colliders 
is being done by comparing a e jf p to the derived limits on a™ m and a™ ax as described 
above. Generally speaking, there are three possible outcomes of this test, which we 
shall now consider in turn. 

1. The test will be deemed successful if it turns out that 

vain , exp -max In 1\ 

°A <°A A ■ K 2 - 1 ) 

Under those circumstances, the dark matter particle discovered at colliders can 
solely account for all of the dark matter in the universe, i.e., there is no need 
for another independent dark matter candidate. This is the scenario which 
has attracted the most attention in the literature. The gap between <r™ n and 
a max j g i nc [i ca tive of the precision with which colliders can test the WIMP 
dark matter hypothesis, and with the availability of a next generation lepton 
collider, typically the amount of wiggle room between a™ m and a™ ax is reduced 
to the order of the current uncertainty in a e ^ p . If this turns out to be the case, 
we will have a triumphant confirmation of the WIMP dark matter hypothesis, 
although one could always find some caveats. 1 

Our major point here is that even though this case has been most widely 
considered in the literature, one should still pay proper attention to the other 

1 For example, it would be very interesting to see how one could reconcile such a successful test 
of dark matter at colliders with a possible experimental direct detection of dark matter in a mass 
range different from the mass range found at colliders. 
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two possible outcomes below, especially since they would clearly indicate the 
presence of some additional new physics. 



2. It is also quite possible that the collider test of dark matter will reveal that 



Under the conventional interpretation, the dark matter candidate would then 
overdose the universe and cannot constitute the present-day cosmological dark 
matter. This conclusion is also subject to caveats either on the astrophysics 
or particle physics side. First, nonstandard cosmological evolution will change 
the picture - for example late inflation may dilute the dark matter density so 
that the particle observed at colliders is indeed the dark matter. This scenario 
can be strengthened by the observation of a direct detection signal of dark 
matter with mass in the range observed at colliders. The information obtained 
at colliders will then provide invaluable insights into early universe cosmology. 
Alternatively, one may look for particle physics resolutions of the puzzle. One 
possible explanation is that the WIMP dark matter candidate found at col- 
liders is metastable and decays post-freeze out to a lighter particle with only 
gravitational interactions (superWIMP), which is the real dark matter [54,55]. 
The mass ratio of the WIMP and the superWIMP can be adjusted to com- 
pensate for the overclosure of the universe and obtain the proper relic density 
today, although significant constraints arise from BBN and large scale struc- 
tures. Colliders are very useful in providing information, such as the mass of 
the dark matter superWIMP, which cannot be obtained by any other means. 
Indeed, direct and indirect detection experiments are bound to fail in their 
searches for purely gravitationally interacting particle dark matter. Another 
possible particle physics explanation is to invalidate the bound a™ ax < a e ^ p . 
For example, additional invisible annihilation channels not easily revealed at 
colliders would push the bound on a™ ax higher. In summary, outcome ( |2.2j ) 
requires either a reconcilation of the collider results by modifying standard 
cosmology and/or postulating new physics, or an invalidation of the colllider 
results by postulating new physics. New physics is thus expected either on the 
cosmology or particle physics side. 

3. The third possible outcome of the dark matter test at colliders is that 



The conventional interpretation of this outcome is that the dark matter parti- 
cle observed at colliders is not the only component of the dark matter in the 
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universe, and one must look for another, yet unobserved, dark matter parti- 
cle. Here again one may look for ways to circumvent this conclusion. On the 
astrophysics side, a nonstandard cosmological history which leads to an en- 
hancement of the relic abundance of the dark matter candidate, as is provided 
by the kination-dominated scenarios considered in this paper, can naturally 
accommodate this apparently unsuccessful dark matter particle candidate. On 
the particle physics side, one would have to invent another, independent dark 
matter candidate, which would provide the missing component of the dark mat- 
ter. While this is possible in principle, it requires nonminimal model building - 
the presence of multiple dark matter candidates would require the existence of 
multiple conserved symmetries to ensure separately the stability of these parti- 
cles on cosmological scales. Just as before, one may think of ways to invalidate 
the problematic bound a^ p < <r™ m . However, and this is the major difference 
between outcomes ( |2.2| ) and ( |2.3| ) from a particle physicist's perspective, it is 
clearly almost impossible to invalidate the lower bound on the annihilation 
cross-section, since, as elaborated above, it is derived from processes which 
have already been firmly observed at colliders, so any channels which may 
have been missed, would only increase <r™ n and thus make the discrepancy 
worse. The outcome ( |2.3|) is also very exciting for an astrophysicist, since, in 
the absence of alternative particle dark matter candidates, it would provide 
a direct indication of a cosmological relic abundance enhancement mechanism 
(such as kination domination) and give reasons to hunt for further correlated 
astrophysical/cosmological signatures [44,56]. 

The above discussion underscores the potential importance of high energy collid- 
ers for the understanding of our universe. Colliders have traditionally been viewed 
as tools for proving that the WIMP particle indeed constitutes the dark matter. Now 
we also see that a potential unexpected outcome of the collider tests of the WIMP 
dark matter hypothesis may provide precious insights into early universe cosmology. 

At the same time, one should not undervalue the potential significance of experi- 
ments dedicated to direct dark matter searches. A positive signal in any one of them 
would have several important implications. First and foremost, it would mean that 
the dark matter is real and would eliminate simple alternative explanations such as 
MOND or modified gravity. Of course, modified gravity can still play a role in the 
dark matter story in the second and third possible outcomes of the dark matter test 
at colliders stated above (see e.g. [45,57]), but particle dark matter must play a sig- 
nificant role in gravitational clustering and galaxy formation. Second, it would boost 
the WIMP dark matter hypothesis, as it can (depending on the inferred relic density) 
rule out particles with purely gravitational interactions as a significant component of 
the dark matter halo (e.g. [54,55,58,59]). Finally, as emphasized in [10,18,19], direct 
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detection experiments will provide an important piece of the dark matter puzzle - 
an independent estimate of the mass of the dark matter particle, which can then be 
contrasted with the analogous measurement at colliders. In the spirit of our earlier 
discussion, here again we can consider three possible outcomes of this comparison: 

1. It may turn out that the WIMP candidate found at colliders is lighter than 
the directly detected dark matter particle. This scenario would lead to several 
interesting possible interpretations. For example, the assumption that the dark 
matter is made up of a single component would imply that the particle found 
at colliders cannot be a dark matter particle and would have to decay outside 
the detector. Alternatively, the dark matter may be made up of several com- 
ponents, allowing for the coexistence of several particles as stable dark matter 
candidates. The collider particle may even be the dominant component if its 
direct detection rates are suppressed. In any case, the direct detection results 
would then provide a rough target for the next energy scale which would need 
to be reached and probed by collider experiments. 

2. Conversely, it may turn out that the WIMP candidate found at colliders is 
heavier than the directly detected dark matter particle. Within the single 
component dark matter scenario, the most natural explanation of the discrep- 
ancy would be that the collider WIMP decays invisibly to the dark matter 
particle. Then, it would be interesting to go back and scour the collider data 
for events where the lighter dark matter particle can directly manifest itself. 

3. Finally, the collider and astroparticle mass determinations may turn out to be 
in agreement, which would point towards a single WIMP dark matter compo- 
nent. The collider measurements of the dark matter particle properties can 
then be used to reconstruct the WIMP annhilation rate in the early universe. 
The outcome of this exercise is extremely interesting. For example, either 
a max < a exp Qr ^exp ^ a min wou } c [ necessarily require nonstandard cosmology. 
If a™ ax < (T^ p , the universe appears overdosed and some mechanism of late 
entropy production (such as a phase transition or a late particle decay) is re- 
quired. If, on the other hand, a^ p < a™ m , the dark matter abundance would 
require some kind of a boost, such as a period of kination domination as will 
be discussed in this paper. 

To summarize, the outcome of experimental tests of the WIMP hypothesis will 
provide significant insights into new physics. Although most studies of this issue 
focus on scenarios in which standard cosmology holds and the thermal WIMP dis- 
covered at colliders is the dominant component of the cold dark matter, such that 
a min < a ex P < a max ^ alternative outcomes also warrant careful consideration. Such 
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alternative outcomes can have important implications for new particle physics, such 
as the possibility of several distinct constituents of the cold dark matter, and/or cos- 
mology, such as the possiblity of cosmological enhancement or dilution mechanisms 
for the Hubble expansion of the universe. We focus here on a particular cosmological 
enhancement mechanism for the dark matter relic abundance due to the presence 
of dark energy in the form of quintessence. In this situation, a e A xp < a™ m < a™ ax , 
and yet the thermal WIMP measured at colliders can be the dominant component 
of the cold dark matter. Kination-dominated quintessence generically can be tested 
through a number of correlated cosmological implications and signatures (see [56] for 
a discussion and further details), but is also a framework in which dark energy can 
make itself manifest at high energy colliders, in striking contrast to the vast majority 
of known dark energy models. We will now concentrate on this scenario and explore 
the prospects for experimental tests at the LHC and ILC. 

3. Cosmology and the dark matter/dark energy connection 

Given the motivation for considering quintessence models of dark energy which in- 
clude a period of kination domination, we turn to its effects on the dark matter relic 
abundance, as first discussed in [39]. In this section, we provide a self-contained set 
of analytic equations which can be used to map any ordinary dark matter freeze out 
scenario to the dark matter freeze out scenario of kination domination. 2 

In the usual thermal WIMP dark matter scenario (see e.g. [60]), the dark matter 
X is assumed to be initially in chemical equilibrium with the thermal plasma that 
is in equilibrium with the photons. The WIMPs remain in equilibrium as long as 
the annihilation reaction rate = ti x {(Jav) (<ta is the cross section for reactions 
that change the number density n x and the averaging is with respect to a thermal 
ensemble) is much larger than the expansion rate of the universe H = a/ a (a is the 
scale factor). However, when Ta falls below H, n x no longer tracks the equilibrium 
density and n x a 3 is nearly an adiabatic constant. This transition from equilibrium 
abundance tracking to a nearly adiabatic constant behavior is known as the freeze 
out transition. Typically, the energy density that governs H during the freeze out 
period is composed of relativistic degrees of freedom in equilibrium with the photons. 

The assumption that the dark energy is in the form of quintessence naturally 
leads one to question whether the relativistic degrees of freedom governing the ex- 
pansion rate H need be in chemical equilibrium with the photons during WIMP 
freeze out. Indeed, since the coupling of the quintessence field to ordinary matter 

2 Here we will restrict our attention to a flat FRW universe with ds 2 — dt 2 — a 2 (t)dx 2 governed 
by the standard Einstein-Hilbert gravitational action. Planck's constant is defined to be M p i 
1.22 x 10 19 GeV. 
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needs to be very weak to maintain a long dynamical time scale and evade fifth force 
constraints, quintessence naturally never reaches chemical equilibrium with the pho- 
tons. However, if the quintessence field were to dominate the universe during freeze 
out, its energy density must dilute faster than the energy density of the SM relativis- 
tic degrees of freedom pn to evade the BBN constraint on the number of relativistic 
degrees of freedom. This constraint can be expressed (at 95% confidence level) as [61] 

A<7,s<1.5^, (3.1) 

where Ag^s is the number of effective degrees of freedom [60] other than the photons, 
electrons, and neutrinos entering the total entropy density 

2tt 2 

S = ^(10.75 + A^ 5 )T 3 , (3.2) 

where To is the temperature during BBN. 

As pointed out in [39] , one way to achieve this is to have a period of kination 
domination, in which the kinetic energy of the quintessence field <E> dominates the 
energy density: p$ ~ |$ 2 Pr. In this case, the energy density scales as 

— oc — . 3.3 

Pr a 2 

This indicates that if the ratio of energy density of $ and photons at the temperature 
of 1 MeV is 

?7$ = |T=lMeV, (3.4) 

the ratio at the time of freeze out is of the order 

— |r=iGeV ~ 10V (3.5) 
Pj 

This indicates that the $ energy density dominates over that of the relativistic de- 
grees of freedom, since p$ ~ 10 5 i]^p R at the approximate period of freeze out and 

< 77* < 1. (3.6) 

Therefore, we can easily arrange p$ to control the freeze out temperature T K in the 
kination dominated scenario. (Note that one should distinguish between p 7 which is 
the density of photons and pr which is the energy density of all relativistic particles.) 

Following the computational approach of [39], we can easily give a map between 
the dark matter abundance of the standard scenario and that of the kination domi- 
nation scenario. With LSP dark matter in mind, the thermal averaged annihilation 
cross section (defined to be proportional to the number changing reaction rate) is 

(ctav) = a + bx, (3.7) 
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where 



x 



T 



(3.8) 



m x is the mass of the WIMP, and a (b) is the s-wave (p-wave) contribution to the 
annihilation cross section. The freeze out temperature parameter will be denoted 
by xjj in the standard scenario and xk in the kination domination scenario. These 
quantities are defined by the equations 



(a + bx UtK )n eq (x UtK ) = H U)K (x u>K ), 

Xu,K 

n^x)=2ml^e-K 



Hu(x) 



87r3 , \ m x 2 

g*{x)—r-x , 



90 



M, 



P i 



H K (x) 



8tt 3 



m 



x x 2 



90 J ~ K ~ M p i 

where g*(x) is defined by the equation for the radiation energy density 

r 2 



7T 



PR = ^9*{x)x A m 4 x . 



(3.9) 
(3.10) 
(3.11) 

(3.12) 
(3.13) 



In the kination domination scenario, the relic abundance today VLk can then be 
expressed relative to the usual relic abundance Vtjj as 



K _ Hu xfj g*{x V/ 



1 + ax 2 



K 



a + bxjj 
a + bxx 



Pu = 1 + 



a + bxu/2 
xu V a + bxu 



2 

x K 



UK — 1 + — \l 1 + OLX 2 K 



aA{u) + bx K B(u) 
a + bxK 



A(u) = - \n[u + VTTv 2 }, 



u 



B{u) 



VI + u 2 - 1 



U = yfciXK 



a = 77$ 



MT )/2]T 2 



9*{Tk) 



(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 
(3.19) 

(3.20) 
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In the above formulas, a subscript U denotes quantities associated with the usual 
freeze out scenario, while a subscript K denotes quantities associated with the ki- 
nation domination scenario. Here the number of relativistic degrees of freedom 
governing the energy density g* obeys the approximate relationship g* « g„ s (e.g. 
g*(Tx) ~ 90 and g*(T ) « 10.75), and T — 1 MeV is the temperature relevant for 
BBN (as opposed to the temperature today). 

Assuming that r/q, < 1 is sufficiently large such that kination dominates at freeze 
out, let us approximate Eq. ( |3.14|) in the limit of either s-wave (a » bxx) or p-wave 
(a << bxx) dominance to obtain intuition for the type of enhancement obtained 
for kination domination. This is useful because in any given realistic model the 
dark matter candidate tends to be either a predominantly p-wave annihilator (e.g. 
the neutralino in supersymmetry) or an s-wave annihilator (e.g. the dark matter 
candidate in most other cases). For a p-wave annihilator we obtain 

V K _.9*s(Tu) T 2 (321) 



Ma g* s (T ) T K T yj2g«{Tu) 

The powers of the temperature T in Eq. ( |3.21| ) can be understood as follows. Since 
&u,k oc (fl(7,x) 3 /atoday an d the volume dilution behaves as 1/a 3 oc T 3 , the volume 
factors alone contribute (Tjj/Tk) 3 - Residual annihilations after freeze out contribute 
a factor of order Tk/Tjj for p-wave dominance. The most nontrivial aspect of the 
kination scenario is that since the freeze out condition is (av)n x = m x H/T, with 

BK-yft (3-22) 



T Q J M pl 

(recall H 2 oc r?$/a 6 during kination domination), Hjj ~ Ty/M p i, and (av) oc T {p- 
wave dominance), there is a factor of Tk/Tq. Note that since Tq/Tu ~ 10" 3 and the 
difference between Tk and Tjj is only logarthmically dependent on the $ 2 energy 



density, it is unrealistic to obtain a ratio smaller than unity in Eq. ( 3.21 



We similarly approximate Eq. (|3.14j) in the limit of s-wave dominance to obtain 

~ g*s(T ) T Q \ og (2u)^gJJ\j)/2 

Inserting typical numbers {m x ~ 100 GeV, T ~ 10~ 3 GeV, g*s(Tu) ~ 9*{Tk) ~ 
10 2 , g» s (T ) « 10.75, b = 5.4 x 10~ 8 GeV~ 2 ,7/$ ~ 1} into Eq. (ggg) and Eq. Qpg ), 
we generically find for both cases 

^ ~ 10 3 . (3.24) 

Hence, the kination scenario relic abundance is much larger than the usual freeze 
out scenario relic abundance for the same microphysical parameters governing the 
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Figure 1: The dark matter relic abundance Qxh 2 in the kination domination scenario, 



as a function of the kination parameter (3.4), for the case of a pure s-wave annihilator 
(d = 0.8 pb, b = 0) or a pure p-wave annihilator (d = 0, b = 36 pb) with spin 1/2 and 



mass m x = 100 GeV. The (blue) solid lines give the exact result from Eq. ( 3.14\ ), while 
the (red) dashed lines correspond to the approximations ( 3.2c ) and ( 13.21 ), The horizontal 
dotted line and the (yellow) shaded band denote the current central value and 2a range 
for the experimental determination of the dark matter relic abundance. 



WIMP annihilation cross section. This is illustrated in Figure ffl, which shows the 
dark matter relic abundance Vl^h 2 in the kination domination scenario, as a function 
of the kination parameter ( |3.4|) , for the case of a pure s-wave annihilator (a = 0.8 
pb, b = 0) or a pure p-wave annihilator (a = 0, b = 36 pb). We assume that 
the dark matter particle is a fermion with spin 1/2 and mass m x = 100 GeV. The 
values for its annihilation cross-sections are chosen so that in the usual scenario 
(77$ = 0) this particle would make up all of the dark matter, in agreement with 
the experimental determination of the dark matter relic abundance by WMAP and 
SDSS (Q x h 2 = 0.111^015 at 2a [1], as indicated by the horizontal dotted line and the 
(yellow) shaded band). The (blue) solid lines give the exact result from Eq. ( |3.14j ), 
while the (red) dashed lines correspond to the approximations ( |3.23| ) and ( |3.21| ). We 
see that the approximate expressions work quite well, down to 77$ ~ 10~ 7 for s- 
wave annihilators, and 77$ ~ 10~ 5 for p-wave annihilators. While the approximations 
( |3.23| ) and ( |3.21|) are useful for understanding the scaling of the relic density in the 
presence of quintessence, in our numerical results below we shall always make use 
of the exact expression ( p.l4|) . From Fig. [I] we also see that, in agreement with the 
naive expectation ( |3.24| ), for 77$ ~ 1 and at this value of m x = 100 GeV, kination 
domination provides an enhancement by about three orders of magnitude of the dark 
matter thermal relic density. The kination enhancement of the dark matter relic 
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Figure 2: The dark matter relic abundance Qxh 2 in the kination domination scenario, as a 
function of the mass m x of the dark matter particle, for the case of a pure p-wave annihilator 
(a = 0, b = 36 ph) and different values of the kination parameter 77$. The horizontal 
(yellow) shaded band denotes the current 2a range for the experimental determination of 
the dark matter relic abundance. 

density can be even more pronounced, if the mass m x of the dark matter particle is 
larger, as evidenced in Figure |2|. There we plot the dark matter relic abundance Q^h 2 
in the kination domination scenario, this time as a function of the mass m x of the 
dark matter particle, for the case of a pure p-wave annihilator (a = 0, b = 36 pb) and 
different values of the kination parameter 77$. As in Fig. [1|, the horizontal (yellow) 
shaded band denotes the current 2a range for the experimental determination of the 
dark matter relic abundance. We see for example, that if the dark matter matter 
particle had a mass of order 1 TeV, kination may provide an enhancement of its relic 
density of up to 4 orders of magnitude! 

Given these results, let us return to the issue of testing the WIMP hypothesis at 
colliders as discussed in Section ^. To make the discussion more concrete, consider 
the following scenario. Suppose that colliders have found a WIMP of mass 100 GeV 
whose spin was measured to be 1/2. The lightest neutralino in supersymmetry is a 
standard such example. Furthermore, because of its Majorana nature, the neutralino 
typically annihilates predominantly in a p-wave, so a << b. The measurements of 
the neutralino couplings and the superpartner mass spectrum can then be translated 
into a bound on the annihilation cross-section <ta, in this case its b component. 
Given the bounds on b, one can then reconstruct the allowed range for the WIMP 
relic abundance, both with and without the effects of kination domination, as shown 
in Fig. [| There we plot the dark matter relic abundance Qxh 2 in the kination 
domination scenario, as a function of the annihilation cross-section b of the dark 
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Figure 3: The dark matter relic abundance Qxh 2 in the kination domination scenario, 
as a function of the annihilation cross-section b of the dark matter particle, for the case 
of a pure p-wave annihilator (d = 0) of spin 1 /2 and mass m x = 100 GeV, for different 
values of the kination parameter . The solid lines are plotted for 77$ values as labelled on 
the plot (every other decade), while the dotted lines correspond to intermediate decades. 
The horizontal (yellow) shaded band denotes the current 2a range for the experimental 
determination of the dark matter relic abundance. 



matter particle, for different values of the kination parameter 77$. As in Figs. [I] 
and 0, the horizontal (yellow) shaded band denotes the current 2<r range for the 
experimental determination of the dark matter relic abundance. The prediction of 
standard cosmology (no kination domination) is given by the rj^ = line. 

As we alluded to earlier, the case which is especially interesting is the one where 
b is determined to be too large, e.g. larger than, say, 100 pb, corresponding to out- 
come ( |2.3| ) in our discussion in the previous section. 3 Such a result can be interpreted 
in one of two ways. First, assuming standard cosmology, one can read off from the 
77$ = line the maximum fraction that such a WIMP can contribute to the dark 
matter budget of the Universe. Alternatively, assuming that the WIMP makes up 



3 Indeed, other outcomes of the collider tests of the WIMP hypothesis, Eqs. ( |2.l| ) and (2.2), 
would strongly disfavor the kination-dominated quintessence dark matter scenario and place an 
upper bound on the quintessence parameter 77$. Of course, quintessence die-hards may still argue 
that there are extra unseen annihilation channels into non-SM states. The corresponding increase 
in a a can be compensated by a nonzero value of 77$. While this is possible in principle, it would 
seem coincidental and rather fine-tuned. 
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100% of the dark matter in the universe and following the horizontal shaded band, 
the lower bound on b would imply a lower bound on 1]$, i.e., a minimum value for 
the kinetic energy contribution of the quintessence field $ to the total energy budget 
in the early universe. We would like to emphasize once again that the lower bound 
on b would tend to be rather robust and difficult to invalidate - it would be derived 
based on some observed channel which would guarantee a minimum value for the 
WIMP annihilation rate. 

4. Numerical Results 

In the previous section we have seen the possible interplay between dark matter and 
dark energy in the early universe. By now it is well appreciated that high-energy 
colliders offer the unique opportunity of creating and studying dark matter in the lab. 
Therefore, if there is a cosmological connection between dark matter and dark energy, 
high-energy colliders can also shed light on the nature of dark energy itself. We now 
analyze specific WIMP scenarios in the context of low energy supersymmetry, for 
which the LSP dark matter candidate is the lightest neutralino, to ascertain to what 
extent forthcoming and future collider experiments can probe dark energy and dark 
matter within this class of quintessence models of dark energy. 

4.1 Choice of benchmark models 

Hence, in the remainder of this section we shall use several study cases to investigate 
the capabilities of the LHC and ILC in determining the relevant dark matter and dark 
energy parameters. It has become customary to perform such studies, using specific 
"study points" in the parameter space of simple models. In particular, in case of su- 
persymmetry, the model of choice has been "mSUGRA" , as it has rather few input 
parameters. Three of them fix the values of the soft supersymmetry breaking pa- 
rameters at the unification scale: the universal scalar mass m , the universal gaugino 
mass Mi/2, and the common trilinear term Aq. The remaining two mSUGRA model 
parameters are tan f3, the ratio of the Higgs vacuum expectation values, and the sign 
of the fi term in the superpotential. While the mSUGRA model is not represen- 
tative of every possible incarnation of supersymmetry, it is nevertheless sufficiently 
general to exhibit four different regions in parameter space with a good dark matter 
candidate. Correspondingly, all of the recent sets of benchmark points suggested in 
the literature [62-64] have zeroed in on those regions, and proposed study points 
where the supersymmetric dark matter candidate (neutralino) makes up all of the 
dark matter in the universe. 

In our case, we would like to include in our analysis the effect of kination, which, 
as we have already seen, tends to enhance the nominal prediction of the dark matter 
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relic density. Therefore, the usual sets of mSUGRA benchmark points are ill suited 
for our purposes, if our goal is to provide a measurement of the kination parameter 
at colliders. Indeed, if the nominal calculation already yields a prediction in exact 
agreement with experiment, kination will then make things worse, and we can at best 
only place an upper limit on 77$. We have therefore chosen to modify the original 
set of LCC benchmark points of [18], so that the nominal calculation would yield 
a value for the relic density which is insufficient to explain all of the dark matter 
in the universe. Barring the existence of another, undiscovered yet dark matter 
candidate, the collider results could then be interpreted as measurements of rjq, and 
would provide a non-trivial link between dark matter and dark energy. The values 
for the original LCC benchmark points (LCC1-LCC4) and our modified versions 
(LCC1 / -LCC4 / ) are listed in Table [1|. The motivation behind each choice will be 
discussed in the following subsections. At this point we would only mention that 
we have tried to only minimally deviate from the original LCC points, and that in 
all but one case our point differs from its LCC counterpart in the value of a single 
mSUGRA parameter. 

Previous studies have already estimated the expected accuracy with which high 
energy physics experiments can pinpoint the supersymmetry parameters at the origi- 
nal LCC benchmark points. One should keep in mind, however, that those estimates 
are only "best guesses" so far, and one would have a better idea of the actual preci- 
sion only after the colliders have been operational for some time, which would allow 
for better understanding of the systematic uncertainties. Nevertheless, we can reli- 
ably estimate the collider precision for our LCC benchmarks. The LCC benchmark 
points yield rather large values for the WIMP annihilation rate, which is typically 
due to the dominance of a single channel: either an s-channel resonance or a coan- 
nihilation channel. Under those circumstances, it is important to know with great 
precision only the masses of the particles involved in the dominant channel, while the 
rest of the spectrum may remain rather uncertain, as long as it does not contribute 
significantly to the WIMP annihilation rate. In what follows, we shall therefore base 
our prediction for the collider precision in measuring the cosmological parameters, on 
the expected precision in measuring the masses of the particles entering the dominant 
annihilation channel. A more sophisticated analysis using the full information about 
the SUSY spectrum and utilizing Markov chain probabilistic techniques along the 
lines of [18] is beyond the scope of our paper, and will not, we believe, significantly 
change our conclusions. 

4.2 LCC1': a study point in the bulk region 

In the remainder of this section we shall present our results for the expected accu- 
racy in determining the dark matter properties at colliders. We shall first discuss 
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Table 1: mSUGRA parameter sets for the four LCC study points and their variations used 
in this study. We also show the relic density and the superpartner spectrum (the masses 
are listed in units of GeV) predicted at each point as calculated with micrOMEGAs [65-67] 
and ISA JET version 7.69 [68]. 



the physics and the relevant measurements at each one of our benchmark points 
LCC1'-LCC4'. Our predictions for the expected allowed range in Vt K h 2 and 77$ as 
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determined by the LHC and ILC for each study point will be summarized in Fig. |9[ 
The superpartner mass spectra for the four original LCC study points and for our 
modified LCC study points are listed in Table [I]. 

The LCC1 study point was chosen in the bulk region, where the sleptons are rel- 
atively light, and neutralino annihilation proceeds predominantly through t-channel 
right-handed slepton exchange. At LCC1, there are 7 channels which contribute 
more than 1% to the annihilation rate, but the three dominant channels are t + t~ 
(32%) and and e + e~ (at 29% each). The remaining channels (neutrino pairs 

and bb) are about 2% each. As one can see, the right-handed sleptons play the 
most important role in mediating neutralino annihilations - first, because they are 
the lightest sfermions in the spectrum, and second, because they have the largest 
hypercharge (the LSP is 97.4% Bino at this point and its couplings are proportional 
to hypercharge). Under those circumstances, in order to pinpoint the WIMP relic 
density, one has to measure precisely the masses of the right handed sleptons and 
the LSP, at the same time making sure that the remaining sparticles in the spec- 
trum are relatively heavy. At the LHC, the slepton masses are difficult to measure 
in direct slepton production [69,70], due to the relatively small slepton production 
cross-section and large SM backgrounds. Fortunately, point LCC1 has relatively 
light colored superpartners (squarks and gluino) which can be produced abundantly 
at the LHC and in their cascade decays may yield sleptons indirectly. In particular, 
the decays of x\ are predominantly to right-handed sleptons (since they are the only 
sfermions lighter than x°). The measurements are by no means trivial, since one 
would have to extract the slepton and neutralino masses from a sufficient set of ob- 
served kinematic endpoints. However, the expectations are that one could measure 
the LSP mass to within 5% and the neutralino-slepton mass differences to within a 
few GeV (for more details, see [18]). Of course, since the sleptons at point LCC1 are 
within the kinematic reach of the ILC, their masses (as well as the LSP mass) can 
be determined at the per mil level at the ILC. 

As we discussed earlier, our strategy will be to modify the original LCC bench- 
mark points so that the resulting relic density is too low and can accommodate the 
kination dominated quintessence scenario. Notice that the relic density at point 
LCC1 is already too large 4 (almost double the WMAP value), which appears to 
make our job rather difficult. There are several ways to reduce the value of Qh 2 . For 
example, one could attempt to reduce the slepton masses even further. However, 
this will not help in our case - first, it will eventually take us into the coannihilation 
region discussed later on in Section [4.4j . Furthermore, the value of m,Q at point LCC1 
is already small enough so that the slepton masses are dominated by the radiative 

4 To some extent this is related to ensuring that the Higgs boson mass is above the LEP limit 
(we shall return to this point later on). 
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corrections due to the gaugino terms in the RGE's. The slepton masses are there- 
fore much more sensitive to the parameter Mxii- This is why we choose to modify 
the LCC1 point by changing the value of Mi/ 2 . The resulting variation in the relic 
density Vtjjh 2 as calculated within standard cosmology, is shown in Fig. ^a, for fixed 
m = 100 GeV, A = -100 GeV, tan/3 = 10, n > and m t = 175 GeV. The 
horizontal (green) shaded band denotes the current 2a range for the experimental 
determination of the dark matter relic abundance. The vertical (yellow) shaded band 
on the left is ruled out by the negative chargino searches at LEP. Since the parameter 
Mi/2 controls the values of the gaugino masses, when it becomes too small, there will 
be a light (wino-like) chargino in the spectrum. It should be kept in mind that the 
LEP bound on the chargino mass is also a function of the electron sneutrino mass 
u e : the bound is diluted in the presence of a light sneutrino due to destructive inter- 
ference between the Z/7 and u e mediated diagrams of chargino production. When 
moving in the opposite direction - increasing M1/2 - we encounter a region where the 
lightest slepton (fi) becomes increasingly lighter relative to the lightest neutralino 
(Xi) and eventually becomes the LSP in the region denoted by the vertical (yellow) 
shaded band on the right. 

Fig. |]a exhibits three regions where the relic abundance is reduced below the 
WMAP level. On the right-hand side, at large M1/2 and near the stau LSP limit, 
this is due to stau coannihilations (see Sec. |4.4j) . On the left-hand side, at low Mi/ 2 , 
we see two dips in Vtjjh 2 due to resonant neutralino annihilations through a Z (the 
left dip) and the light CP-even Higgs boson h (the right dip). The Z resonance is 
already inside the excluded region but the h resonance is still allowed. That is where 
we chose our modified LCC1' point, marked by the vertical line labelled LCC1'. 

The SUSY spectrum of point LCC1' is given in Table [|. It is somewhat lighter 
than the spectrum at LCC1, but the hierarchy of states is similar, which allows us to 
assume the same precision in the sparticle mass determinations as for point LCC1. 
The chargino mass is 97.5 GeV, but (as discussed above) is still allowed due to the 
presence of a light electron sneutrino, which degrades the LEP chargino bound. One 
might be worried that the light CP-even Higgs boson mass is also lighter - 108.7 
GeV, which is below the LEP Higgs mass limit. Indeed such a light Higgs boson is 
ruled out. However, the exact value of the Higgs boson mass is not essential for our 
analysis, since we are only interested in the precision with which the masses can be 
determined, rather than their actual values. At the original point LCC1 the problem 
was avoided by considering a heavier spectrum, which resulted in an unacceptable 
value for the relic density. Both of these problems can be simply solved by relaxing 
some of the mSUGRA assumptions, for example, scalar mass nonuniversality in the 
third generation [71] can easily lift the Higgs mass, while gaugino non-universality 
can remove any remaining tension with the light chargino bound from LEP. 
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Figure 4: a) The dark matter relic abundance Q.jjh 2 in standard cosmology, as a function 
of the universal gaugino mass parameter M^ji in minimal super gravity, for fixed rriQ = 100 
GeV, A = -100 GeV, tan/3 = 10, fj, > and m t = 175 GeV. The horizontal (green) 
shaded band denotes the current 2a range for the experimental determination of the dark 
matter relic abundance. The vertical (yellow) shaded band on the left is ruled out by the 
negative chargino searches at LEP. The other vertical (yellow) shaded band on the right 
is ruled out because the lightest superpartner (stau) is charged. The vertical line marked 
LCC1 (LCC1') denotes the M 1/2 value for the LCC1 (LCC1') study point (see Table |). 
b) The same as a), but plotted versus the mass difference 2m^o — m^, which indicates 
the proximity to the light Higgs pole. The vertical light blue (yellow) band indicates the 
expected experimental precision in determining the value of the combination 2m^o — 
at the LHC (ILC). 

Given the SUSY spectrum of point LCC1', we now turn our attention to the 
corresponding prediction for the WIMP relic density at colliders. The most important 
feature of our spectrum is that 2m^o m m^, which allows the neutralinos to annihilate 
very efficiently on the Higgs resonance. Unlike point LCC1, we now find that the 
dominant neutralino annihilation channels are bb (90%) and t + t~ (9%), as expected 
on the h resonance. This implies that the Higgs resonant diagram by far overwhelms 
all other neutralino annihilation processes combined. This can already be guessed 
from Fig. ^a, where the turn on of the resonant h diagram causes a sharp drop in 
the relic density of almost 2 orders of magnitude. The same effect is even more 
evident in Fig. |]b, where we plot our results from Fig. |]a versus the mass difference 
2m^o — rrih, which is an indicator of the proximity to the light Higgs pole. Figs. |]a 
and |]b make it clear that in order to predict the neutralino relic abundance, we 
only need to concentrate on the precision with which the masses m^o and can be 
determined at colliders. At the LHC, nth can be precisely measured in the diphoton 
channel, therefore the dominant uncertainty is still from the determination of m^o, 
which we take to be 5%, as for point LCC1. At the ILC, we take the uncertainty on 
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both m^o and rrih to be 0.05 GeV [18]. The resulting uncertainty on the combination 
2m^o — rrih at the LHC (ILC) is shown in Fig. |]b by the light blue (yellow) band. 
The remaining mass spectrum measurements can be assumed to be similar to those 
for point LCC1, but as far as the relic density is concerned, they become largely 
irrelevant, their only significance being to show that we are sufficiently far away 
from other special regions, e.g. the stau coannihilation region. 

Fig. |]b can be used to directly translate the uncertainty on 2m^o — rrih into the 
corresponding uncertainty on Quh 2 , as calculated in standard cosmology. Turning 
on the effect of kination dominance, we then find the result in Fig. [|a, which shows 
the expected precision in the simultaneous determination of the dark matter relic 
abundance Qxh 2 and quintessence parameter 77$ at the LHC (blue band) and the 
ILC (red band). Fig. [|a illustrates our main point - that in kination dominated 
quintessence scenarios, colliders constrain a combination of the dark matter and dark 
energy properties. In the case of point LCC1', measurements at the LHC alone will 
still be consistent with a WIMP hypothesis within standard cosmology. One would 
really need the precision of the ILC in order to see that within standard cosmology, 
this particular WIMP is not enough to explain all of the dark matter. Alternatively, 
if the WIMP does make all of the dark matter in the universe, the ILC, combined 
with precision cosmology, provide a measurement of the quintessence parameter 77$ 
at the percent level. 

4.3 LCC2': a study point in the focus point region 

In this and the next two subsections, we shall repeat the analysis we have just done 
for the LCC1' bulk point, for our other three study points. Study points in the bulk 
region of mSUGRA are relatively "collider-friendly", in the sense that a relatively 
large number of particles can be produced and studied, and as a result, a multitude of 
measurements can be made. The remaining three study points are chosen in regions 
where typically a smaller set of states is accessible at colliders. Nevertheless, one can 
still obtain a similar precision on the relic density determination, since the limited 
knowledge of the spectrum is partially offset by the fact that only a few channels 
dominate the annihilation rate, and often it is only a specific feature of the spectrum 
which governs the annihilation rate and thus needs to be known precisely. 

The LCC2 study point was chosen in the so called "focus point" region of 
mSUGRA [72,73]. The region is characterized by relatively heavy scalars, which 
alleviates the tension with a number of phenomenological constraints: e.g. flavor 
problem, CP-violation, light Higgs boson mass, proton decay etc. [74-76]. The LSP 
in the focus point region is still predominantly Bino, but has a non-negligible Hig- 
gsino component, which opens up new annihilation channels into gauge and/or Higgs 
bosons [77]. For example, at LCC2 the LSP is 68% Bino, 28% Higgsino, and 4% 
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Figure 5: a) The dark matter relic abundance fti/h? in standard cosmology, as a function 
of the universal scalar mass parameter rriQ in minimal super gravity, for fixed Mi/% = 300 
GeV, Aq = GeV, tan/3 = 10, [i > and m* = 175 GeV. The horizontal (green) shaded 
band denotes the current 2a range for the experimental determination of the dark matter 
relic abundance. The vertical (yellow) shaded band on the right is ruled out from the 
negative chargino searches at LEP. The vertical line marked LCC2 (LCC2 1 ) denotes the 
mo value for the LCC2 (LCC2 1 ) study point (see Table b) The masses of the charginos 
(dashed, red) and neutralinos (solid, blue), and the value of the \x parameter (dotted, green) 
as a function of mo, for the same fixed parameters as in a). 

Wino, and the dominant annihilation channels are W + W~ (78%), ZZ (12%) and 
Zh (6%). The prediction of the relic density in this region is perhaps most sensitive 
to the amount of gaugino-higgsino mixing: increasing (decreasing) the Higgsino com- 
ponent of the LSP decreases (increases) the relic abundance. This is illustrated in 
Fig. |5|a, where we show the dark matter relic abundance Qjjh 2 in standard cosmology 
as a function of the universal scalar mass parameter m , for fixed M\m = 300 GeV, 
Aq = GeV, tan /3 = 10, \i > and m t = 175 GeV. As before, the horizontal (green) 
shaded band denotes the current 2a range for the experimental determination of the 
dark matter relic abundance. The value of the Higgsino mass parameter \x is anti- 
correlated with mo, as seen in Fig. [5|b, where we show the masses of the charginos 
(dashed, red) and neutralinos (solid, blue), and the value of the \i parameter (dot- 
ted, green) as a function of mo, for the same fixed parameters as in Fig. |^a. We see 
that increasing the value of mo away from the LCC2 point will lead to a smaller \x, 
larger Higgsino component in the LSP, and correspondingly a lower Quh 2 . This is 
the motivation behind our choice of the new value of m = 3360 GeV for our modi- 
fied study point LCC2', as marked by the vertical line in Figs. ||a and ||b. At point 
LCC2', the LSP is well mixed: 56% Higgsino, 38% Bino and 6% Wino. The domi- 
nant annihilation channels are still W + W~ (88%) and ZZ (5%). If we increase m 
even further, the LSP becomes purely Higgsino-like, and its mass begins to track the 
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decreasing value of /i, eventually closing the W + W~ and ZZ annihilation channels, 
which leads to the local increase in VLjjh 2 observed in Fig. [|a. While this increase 
is in principle sufficient to bring the prediction for VLuh 2 back to the desired level, 
it takes place inside a region of parameter space (shaded in yellow in Fig. where 
the lightest (Higgsino-like) chargino is too light and is ruled out by LEP. 5 Just as in 
Fig. f|a, inside the light chargino region we again find two sharp dips in Qjjh 2 , which 
correspond to resonant annihilations on the h and Z pole, respectively. 

Points LCC2 and LCC2' have similar spectra. The sfermions are very heavy and 
will escape detection at the LHC and ILC Gluino production at the LHC will lead to 
long decay chains yielding jets, leptons and missing energy. B-tagging can improve 
the sensitivity of the gluino search [79] and allow for a gluino mass measurement at 
the LHC [80]. However, in the focus point region, the relic density is determined 
primarily by the properties of the neutralino and chargino sectors, which will be 
mapped out relatively well at the ILC [12,81,82]. In particular, all but the heaviest 
chargino and neutralino states are accessible to the ILC at a center-of-mass energy of 
500 GeV. Given the expected set of measurements available at point LCC2 [81], [18] 
estimates the resulting uncertainty in VLuti 2 as 82% at the LHC and 14% at the ILC. 
Given that the SUSY spectrum at our point LCC2' has only minor quantitative dif- 
ferences from the spectrum at LCC2, we expect a similar precision on Qjjh 2 at point 
LCC2' as well. Using the same uncertainty levels (82% at LHC and 14% at ILC), 
we show in Fig. [|b the resulting expectations for the simultaneous determination 
of the dark matter relic abundance Qxh 2 and quintessence parameter 7/$ at point 
LCC2' for the LHC (blue band) and the ILC (red band). Our result should be taken 
with a grain of salt, since we have not done a dedicated analysis of the spectroscopy 
measurements at point LCC2'. Nevertheless, it seems that already at the LHC one 
would be able to place a lower bound on the quintessence parameter 77$, and limit its 
value to within 2-3 orders of magnitude. The ILC, in turn, will significantly narrow 
down the allowed range for 77$ and within the context of our scenario, provide a 
unique probe of dark energy, unavailable by other means. 

4.4 LCC3': a study point in the stau coannihilation region 

In most of the mSUGRA parameter space, the LSP is sufficiently lighter than the 
rest of the SUSY spectrum, so that the relic density is determined primarily by the 
LSP self-annihilation rate. However, mSUGRA also exhibits special regions where 
the LSP is sufficiently degenerate with another supersymmetric particle, so that both 
are present at freeze out and can affect the resulting dark matter relic density both 
through their self-annihilations as well as their co-annihilations with each other. In 

5 The chargino constraint can be evaded for larger values of mo and M 1 / 2 , where the pure Higgsino 
LSP can become a viable dark matter candidate [78] . 
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general, turning on coannihilations can lead to an increase or a decrease of the relic 
density, depending on the type and properties of the coannihilating particle. In 
mSUGRA the LSP is usually a Bino-like neutralino, whose self-annihilation rates 
are typically rather small, overdosing the universe. The presence of coannihilations 
would then typically tend to enhance the overall effective annihilation rate and lower 
the dark matter relic abundance. 

Point LCC3 was chosen in the region where the lightest neutralino is very close 
in mass to the lightest tau slepton f x (see Table 0). Their mass splitting was carefully 
adjusted so that the neutralino-stau coannihilations [83] would dilute the relic density 
precisely to the WMAP levels. In addition to the usual neutralino annihilation 
channels x\x\ - bb (21%), xlxl - t+t' (12%), - and - e+e" 

(5% each); we also have sizable coannihilation effects: XiTi ~~ > hr (21%), Xi^i ~^ 1 T 
(17%) and Xi T i — * Zt (6%). Naturally, reducing the neutralino-stau mass splitting 
even further would enhance the coannihilation contribution, and drop Qjjh 2 below 
WMAP levels, which is what we need for our modified study point. A simple way to 
control the Xi'^i mass splitting is provided by the parameter m which affects the 
mass of fi, but not Xv This is illustrated in Fig. |]a where we show the dark matter 
relic abundance Qjjh 2 in standard cosmology, as a function of the universal scalar 
mass parameter mo in minimal supergravity, for fixed M\/2 = 360 GeV, Aq = GeV, 
tan p = 40, fi > and m t = 175 GeV. The horizontal (green) shaded band denotes 
the current 2a range for the experimental determination of the dark matter relic 
abundance. The vertical line marked LCC3 (LCC3') denotes the mo value for the 
LCC3 (LCC3') study point (see Table []]). The value of mo cannot be too low, since 
then fi becomes the LSP (in the vertical yellow-shaded band). Nevertheless, close to 
the charged LSP region, the relic density gets to about an order of magnitude below 
the WMAP preferred value. At point LCC3', the neutralino-stau mass splitting is 
reduced from 12.3 GeV down to 4.7 GeV. This is sufficient to make coannihilations 
and stau annihilations dominate the relic density calculation. At point LCC3' we find 
that XiXi annihilation processes only account for about 6% of the total annihilation 
rate. The dominant channels for dark matter number- changing processes are Xi^i ~~ > 
hr (25%), x\h -> ir (17%), ffrf -> hh (17%), r+ff -> t+t~ (11%), etc. 

Precision spectroscopy at points LCC3 and LCC3' is rather challenging at the 
LHC. For one, the leptonic decay chain X2 ^ + ^~Xi is lost, as x\ predominantly 
decays to taus. Furthermore, because of the small mass splitting, the taus 

from the f\ decay tend to be relatively soft and difficult to reconstruct. A couple of 
recent analyses have attempted simultaneous extraction of the neutralino-stau mass 
difference and the gluino mass in a channel with jets, 2r's and missing energy [84] 
or a channel with jets, 3r's plus missing energy [85]. However, it was found that 
those methods fail at very low neutralino-stau mass differences - below 5 GeV, as is 
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Figure 6: a) The dark matter relic abundance Qjjh 2 in standard cosmology, as a function of 
the universal scalar mass parameter mo in minimal super gravity, for fixed M\j2 = 360 GeV, 
Aq = GeV, tan/3 = 40, fi > and m t = 175 GeV. The horizontal (green) shaded band 
denotes the current 2a range for the experimental determination of the dark matter relic 
abundance. The vertical (yellow) shaded band on the left is ruled out because the lightest 
superpartner (stau) is charged. The vertical line marked LCC3 (LCC3') denotes the mo 
value for the LCC3 (LCC3') study point (see Table [J), b) The same as a), but plotted versus 
the mass difference ~ m x1 > which controls the effect of stau coannihilations. The vertical 
light blue (yellow) band indicates the expected experimental precision in determining the 
value of the mass splitting — m^o at the LHC (ILC). 



the case of our point LCC3'. We shall therefore conservatively assume that in the 
absence of a definitive measurement, the LHC can only rule out large enough mass 
splittings (~ 40 GeV) which would have made a measurement possible. The ILC, on 
the other hand, can measure a neutralino-stau mass splitting as low as 5 GeV, down 
to about ±1 GeV [86,87]. These uncertainties are used in Fig. |b, where we plot the 
result from Fig. |6]a versus the mass difference — "^x?' wn i cn controls the effect 
of stau coannihilations. The vertical light blue (yellow) band indicates the expected 
experimental precision in determining the value of the mass splitting — m^o at 
the LHC (ILC), as discussed above. We can now use the result from Fig. ||b to 
anticipate the precision (shown in Fig. ||c) in the simultaneous determination of the 
dark matter relic abundance Qxh 2 and quintessence parameter 77$ at point LCC3', 
for the case of the LHC (blue band) and the ILC (red band). The challenges at 
the LHC mentioned earlier are readily evident, as the LHC results alone will still 
be consistent with standard cosmology, ruling out only the largest possible values of 
77$. With the addition of the ILC, one is again able to pinpoint quite accurately the 
quintessence parameter, and rule out the WIMP scenario within standard cosmology. 
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4.5 LCC4': a study point in the Higgs funnel region 

Our final example is a study point illustrating the so called "Higgs funnel" case 
[88,89]. The SUSY spectrum at point LCC4 is somewhat similar to the one at 
point LCC3, however, the relic density is controlled by a different physics process. 
Here neutralino annihilations occur sufficiently close to the A resonance, so that 
the relevant parameters are the heavy Higgs masses and widths. In this sense, the 



situation is similar to our modified "bulk" point studied in Sec. [L|, where neutralinos 
annihilated near the light Higgs pole. 

When neutralino annihilations occur exactly on the A pole, the relic density is 
typically too low, since the annihilation cross-sections are suppressed only by the 
Higgs width instead of the Higgs mass. Conversely, when neutralinos annihilate far 
away from the resonant pole, the relic density is determined by the other annihilation 
channels and is typically too large (unless we are in one of the regions discussed in 
the previous three sections). Therefore, there are two special places, on both sides of 
the A pole, where the relic density will be just right. The parameters for point LCC4 
were chosen so that neutralino annihilations take place in the right place below the 
A resonance. The dominant annihilation channels at point LCC4 are X1X1 ~^ bb 
(78%) and x\x\ t + t~ (14%), indicative of the heavy Higgs branching fractions 
into fermion pairs. 

Following our usual procedure, we wish to modify the LCC4 point so that to 
reduce the relic abundance. Since we are already in the vicinity of the A pole, 
the simplest way to achieve this without leaving the funnel region, is to simply 
reduce the heavy Higgs masses so that neutralino annihilations proceed on resonance. 
In mSUGRA, the Higgs masses are directly controlled by the mo parameter and 
therefore can be reduced by lowering m . Fig. [j] shows the resulting variation of 
the dark matter relic abundance Vljjh 2 in standard cosmology, for fixed values of the 
remaining LCC4 parameters (M 1 / 2 = 420 GeV, A = GeV, tan/5 = 53, fi > and 
m t = 178 GeV). As usual, the horizontal (green) shaded band denotes the current 
2(7 range for the experimental determination of the dark matter relic abundance. 
The vertical line marked LCC4 denotes the mo value for the LCC4 study point (see 
Table [I]). One might attribute the observed sharp reduction in the relic density at 
lower mo to approaching the heavy Higgs pole, however this is not really the case. 
In the vertical (yellow) shaded band on the left we find a stau LSP, hence this region 
is cosmologically ruled out. Furthermore, near its boundary the neutralino and stau 
are quite degenerate, and we encounter the stau coannihilation situation discussed in 
Sec. [14| . We have checked that the reduction in VLjjh 2 near the stau LSP boundary 
is primarily due to stau coannihilations as opposed to resonant annihilations on the 
Higgs pole (on the stau LSP boundary we still find m^ — 2m^o ~ 75 GeV). Therefore, 
the m variation in the vicinity of the LCC4 study point does not reveal the classic 
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Figure 7: The dark matter relic abundance Quh 2 in standard cosmology, as a function 
of the universal scalar mass parameter mo in minimal super gravity, for fixed M]/ 2 = 420 
GeV, Aq = GeV, tan/3 = 53, [i > and mt = 178 GeV. The horizontal (green) shaded 
band denotes the current 2a range for the experimental determination of the dark matter 
relic abundance. The vertical (yellow) shaded band on the left is ruled out because the 
lightest superpartner (stau) is charged. The vertical line marked LCC4 denotes the mo 
value for the LCC4 study point (see Table [J). 



two-sided funnel shape, since the other half of the funnel is obscured by the stau 
coannihilation region. Similar conclusions hold if we vary tan/3 instead - just like 
mo, tan/? affects in a similar way both the Higgs and stau masses. We therefore 
choose to select our point in a funnel region which is sufficiently far away from the 
stau coannihilation boundary. This can be simply achieved by switching the sign 
of the \x parameter, which reverses the sign of the tan /3 enhanced Yukawa coupling 
corrections, leading to larger stau masses relative to the heavy Higgs masses. We 
therefore modify the values of two parameters (tan/3 = 50 and n, < 0) and in 
Fig. |]a once again we show the variation of the dark matter relic abundance VLjjti 1 
in standard cosmology, as a function of the universal scalar mass parameter mo 
in minimal supergravity. We now clearly observe the classic Higgs funnel shape. 
The minimum of the relic density is found right on resonance, around mo = 940 
GeV, and the WMAP values (shown by the horizontal green shaded band) can be 
achieved on either side. The stau LSP region is now sufficiently far away - in fact 
the vertical yellow-shaded region on the left in Fig. [|a is now ruled out because of 
a light Higgs boson in the spectrum. For our modified study point we choose the 
value m = 950 GeV (denoted by the vertical blue line) where we similarly find the 
dominant annihilation channels to be XiXi ~ * bb (91%) and XiXi ~ > t + t~ (9%). 
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Figure 8: a) The dark matter relic abundance Vtjjh 2 in standard cosmology, as a function 
of the universal scalar mass parameter rriQ in minimal super gravity, for fixed M1/2 = 420 
GeV, Aq = GeV, tan/3 = 50, \x < and mt = 178 GeV. The horizontal (green) shaded 
band denotes the current 2a range for the experimental determination of the dark matter 
relic abundance. The vertical (yellow) shaded band on the left is ruled out because there 
is a very light Higgs boson in the spectrum. The vertical line marked LCC4' denotes 
the mo value for the LCC4' study point (see Table |IJ). b) The same as a), but plotted 
versus the mass combination 2m^o —m^, which indicates the proximity to the heavy Higgs 

-VI 

pole. The vertical light blue (yellow) band indicates the expected experimental precision 
in determining the value of the combination 2m^o — vha at the LHC (ILC-1000). 

We are now in position to discuss our expectations for the precision of the dark 
matter and dark energy determinations at the LHC and ILC at point LCC4'. As 
with any such Higgs funnel point, it is most important to measure the heavy Higgs 
spectrum in addition to the LSP mass. For consistency, we shall again use the 
assumptions of [18], that at point LCC4 the LHC will be able to determine the 
neutralino (heavy Higgs) mass to within 17 GeV (1.5 GeV). Clearly, the error on the 
relevant mass combination 2m%o — itia is then dominated by the error on the LSP 
mass determination. Therefore, we can use the same error estimates for our modified 
point LCC4', which has a very similar value of the LSP mass. In Fig. |8|b we again 
plot the dark matter relic abundance, for the same fixed parameters as in Fig. |8|a, but 
this time versus the relevant mass combination 2m^o — rriA (compare to Fig. |]b). The 
vertical light blue band indicates the expected experimental precision in determining 
the value of the combination 2m^o — rriA at the LHC. Unfortunately, at both point 
LCC4 and LCC4' the heavy Higgs spectrum is too heavy to be observed by the ILC 
at center-of-mass energy of 500 GeV (ILC-500). In addition, at our modified point 
LCC4' the staus are also heavy enough to escape detection at the ILC-500. Since 
ILC-500 does not add anything new to the dark matter determinations, following [18] 
we choose to consider the ILC upgrade at center-of-mass energy of 1000 GeV (ILC- 
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Figure 9: Expected precision in the simultaneous determination of the dark matter relic 
abundance ^Ixh 2 and quintessence parameter rj§ at the LHC (blue band) and the ILC (red 
band), for a) LCC1', b) LCC2I ', c) LCC3' and d) LCC4' study point. 



1000). There the chargino-neutralino sector becomes accessible and the LSP mass 
is expected to be measured to within 1.4 GeV. The vertical yellow band in Fig. ^b 
then indicates the corresponding experimental precision at ILC-1000 in determining 
the value of the combination 2m^o — at our modified point LCC4'. 

The experimental precision in determining the relevant sparticle properties from 
Fig. ^b can now easily be translated into the corresponding uncertainties on the dark 
matter and dark energy properties. The result for our modified point LCC4' is shown 
in Fig. |^d. We see that, as in the case of LCC2', the LHC can already uncover the 
mismatch between the measured WIMP particle properties and those required by 
standard cosmology. The LHC can also provide a surprisingly good determination 
of the kination parameter 77$ - in fact, as can be seen by comparing the thickness of 
the LHC bands in the four panels of Fig. || point LCC4' is where the LHC does the 
best job. The advantage of the ILC as a precision machine is also readily seen in all 4 
panels. For point LCC4' where the ILC-500 energy is not sufficient, it should be kept 
in mind that the ILC project will go ahead only after the first LHC results become 
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available. Once the LHC observes the A resonance near 350 GeV, and indicates that 
all electroweak and colored superpartners are rather heavy, it would be clear that 
the ILC design effort would shift towards the higher energy ILC option. 

5. Discussion and Conclusions 

We are entering a new era in which the discovery of new fundamental physics at 
colliders may help us answer some of the most important puzzles of cosmology. In 
particular, as the LHC is likely to lead to the discovery of WIMPs that are identical 
to the particles which make up dark matter galactic halos, we may soon learn about 
new microphysics which governs the evolution of the universe. Therefore, it is natural 
to ask which aspects of cosmology can be probed by the anticipated identification 
and microphysics of the dark matter particles. 

One answer is that colliders will probe the period of our history at temperatures 
of order of a few GeV when the dark matter particles (assuming they are thermal 
WIMPs) presumably fell out of chemical thermal equilibrium. Therefore, we may be 
able to use collider information to probe properties of every significant energy com- 
ponent in the universe at the WIMP freeze out temperature, including field degrees 
of freedom such as quintessence, which may be responsible for the dark energy today. 
Since quintessence cannot be directly probed at any foreseeable terrestrial controlled 
experiments, we have investigated to what extent LHC and ILC can probe the cos- 
mological properties of a very broad class of quintessence models parameterized by 
a single parameter 1]$ < 1 (the ratio of energy density of quintessence to photon 
energy density at the time of BBN defined as T = 1 MeV), assuming that there was 
a period of kination domination during the time of freeze out. 6 

To compare the power of LHC and ILC in probing the early universe dark energy 
cosmology, we have considered mSUGRA models with parameters analogous to the 
set of four LCC study points [18] which were chosen to represent the four "good" 
dark matter regions. We find that for the cases in which the dark matter annihilates 
primarily through the lightest Higgs resonance and for the stau coannihilation region, 
the LHC is only able to put an upper bound on the parameter 77$ characterizing the 
quintessence dynamics, while the ILC can put both an upper and a lower bound 
on 1]$,. Since a lower bound means a discovery of nontrivial dark energy dynamics 
while an upper bound is consistent with no dark energy component, the ILC has 
the potential to find supporting evidence for the possible kination history of dark 

6 A period of kination domination during which the scalar field energy density dilutes as a~ 6 is 
natural for ?y$ ~ 0{1) because of BBN constraints. The particular parametrization assumes that 
kination domination continues until temperatures of 1 MeV, but the results can be applied even to 
models for which this does not occur by making a model-dependent rescaling of the results. 
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energy. Indeed, it is remarkable that in the coannihilation region, the ILC can 
measure 77$ as small as 1CF 6 , while the LHC can only put an upper bound of 1CF 4 
on the same quantity. This should serve as yet one more important motivation for 
the construction of the ILC. 

Furthermore, such conjectures about the underlying dark energy dynamics have 
other independent observable signatures which can be tested by near-future non- 
collider experiments. For example, observable signatures and correlated constraints 
were discussed in a previous publication [56] for the broad class of kination-dominated 
quintessence scenarios with an embedding in inflationary cosmology relevant for our 
present study. Therefore, just as big-bang nucleosynthesis has served as a sturdy 
pillar of cosmology to constrain dynamics of many conjectures of physics beyond 
the Standard Model, the anticipated identification of dark matter at colliders and 
astrophysics experiments will provide another robust pillar for cosmology, which will 
have profound consequences for uncovering the elusive nature of dark energy. 
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